BUCKLING BEHAVIOR OF COMPRESSION-LOADED QUASI-ISOTROPIC CURVED 

PANELS WITH A CIRCULAR CUTOUT 


Mark W. Hilburger 
National Research Council 
Washington, D. C. 20418 

Vicki O. Britt 

Gulf Stream Aerospace Corp. 
Savannah, GA 31402-2206 

Michael P. Nemeth 
NASA Langley Research Center 
Hampton, Virginia 23681-2199 


Presented at the 40th AIAA/ASME/ASCE/AHS/ASC Structures, 
Structural Dynamics, and Material Conference 

A1AA Paper No. 99-1279 


St. Louis, Missouri 
April 12-15. 1999 



BUCKLING BEHAVIOR OF COMPRESSION-LOADED QUASI-ISOTROPIC CURVED 
PANELS WITH A CIRCULAR CUTOUT 


Mark W. Hilburger* 
National Research Council 
Washington. D.C. 20418 

x 

Vicki O. Britt 

Gulf Stream Aerospace Corp. 
Savannah, GA 31402-2206 

Michael P. Nemeth* 
NASA Langley Research Center 
Hampton, Virginia 23681-0001 


Abstract 

Results from a numerical and experimental study of 
the response of compression-loaded quasi-isotropic 
curved panels with a centrally located circular cutout are 
presented. The numerical results were obtained by using 
a geometrically nonlinear finite element analysis code. 
The effects of cutout size, panel curvature and initial geo- 
metric imperfections on the overall response of compres- 
sion-loaded panels are described. In addition, results are 
presented from a numerical parametric study that indicate 
the effects of elastic circumferential edge restraints on the 
prebuckling and buckling response of a selected panel and 
these numerical results are compared to experimentally 
measured results. These restraints are used to identify the 
effects of circumferential edge restraints that are intro- 
duced by the test fixture that was used in the present study. 
It is shown that circumferential edge restraints can intro- 
duce substantial nonlinear prebuckling deformations into 
shallow compression-loaded curved panels that can re- 
sults in a significant increase in buckling load. 

Introduction 

A common structural element that is found in many 
aerospace structures is the thin-walled curved panel with 
a centrally located circular cutout. Cutouts commonly ap- 
pear in curved panels as access ports, doors, or windows. 
During operation, these structural elements may experi- 
ence compression loads, and thus their buckling and post- 
buckling behavior are important factors in determining 
safe operating conditions and effective designs for these 
structures. 

Several studies have been presented that show that a 
circular cutout in a compression-loaded composite flat 
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plate can have a significant effect on the buckling and 

i •s 

postbuckling response of the structure. However, sig- 
nificantly fewer studies have been presented for com- 
pression-loaded curved panels with a cutout. More 
specifically, limited results on the effects of a circular 
cutout on the buckling response of compression-loaded 
aluminum curved panels are presented in Ref. 4 and lim- 
ited results on the effects of rectangular cutouts and lam- 
inate stacking sequence on the buckling behavior of 
compression-loaded curved composite panels are pre- 
sented in Refs. 5 through 8. In Ref. 4 it has been shown 
that circular cutout size greatly affects the buckling char- 
acteristics of compression-loaded curved isotropic pan- 
els. In particular, the results indicate that aluminum 
panels with cutout-diameter-to-panel-width ratios d/W 
ranging from 0 to 0.3 exhibit nonlinear behavior that in- 
cludes a distinct buckling event with an unstable post- 
buckling response while panels with values of d/W 
ranging from 0.4 to 0.6 exhibit a stable monotonically in- 
creasing nonlinear load-shortening response. These re- 
sults also indicate that as the value of d/W increases, the 
magnitude of the nonlinear prebuckling out-of-plane de- 
formations increases, which results in an increase in the 
nonlinearity of the prebuckling load-shortening response 
curve. 

A review of the results in Refs. 1 through 6 indi- 
cates that the buckling and postbuckling behavior of 
compression-loaded flat plates with a circular cutout is, 
for the most part, well understood. However, the interac- 
tion between circular cutout size, panel curvature, mate- 
rial orthotropy and anisotropy, initial geometric 
imperfections and boundary conditions, and their effects 
on the buckling and postbuckling response of compres- 
sion-loaded composite curved panels are not well under- 
stood. 

Because the amount of information on this impor- 
tant subject is so limited, the results from a numerical and 
experimental investigation of the response of compres- 
sion-loaded | ±45/0/90]^ quasi-isotropic curved pan- 
els with a centrally located circular cutout are presented 
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herein. The objective of the present paper is to identify 
the effects of cutout size, panel curv ature, initial geomet- 
ric imperfections and selected boundary conditions on 
the compression response of these common laminated 
composite structures. Toward that objective, numerical- 
ly predicted and experimentally measured results are 
presented and compared for 15- and 60-in-radius panels 
with cutout-diameter-to-plate-width ratio d/W ranging 
from 0 to 0.6. In addition, numerically predicted results 
are presented that show the effect of circumferential re- 
straint of the loaded edges on the severity of nonlinear 
prebuckling deformations and the buckling load level. 
First, the finite-element models and analysis methods are 
described. Next, the test specimens and test procedure 
are presented. Then, an in-depth discussion of the results 
is given. 

Finite-Element Model and Analysis Methods 

A typical finite-element model of a curved panel 
with a centrally located circular cutout is illustrated in 
Fig. 1. Points on the panel mid-surface are located by an 
x-y-z curvilinear coordinate frame whose origin is at the 
top corner of the panel. The panel length, w idth and ra- 
dius are defined as L. W and R, respectively, and the cut- 
out diameter d is defined as the diameter of a right 
circular cylinder that intersects the panel and whose axis 
is perpendicular to the tangent plane at the center of the 
panel. Panels with 15- and 60-in. radii were studied in 
the present investigation. The panel lengths and arc- 
widths were 14.75 in. and 14.5 in., respectively, for all 
panels. Cutout diameters ranged from 0 to 8.4 in. The 
panels were modeled as 24-ply-thick [±45/0/90]^ 
quasi-isotropic graphite-epoxy laminates with a nominal 
ply thickness of 0.005 in. and a total, nominal panel 
thickness of 0.12 in. The nominal lamina material prop- 
erties are as follows: longitudinal modulus Ej = 17.5 
Msi, transverse modulus E 2 = 1.51 Msi, in-plane shear 
modulus Gj 2 = 0.78 Msi, and major Poisson’s ratio v 12 = 
0.295. Idealizations of the test fixture support conditions 
were used in the finite-element model. To simulate 
clamped loaded edges, the circumferential and out-of- 
plane displacements, v and w, respectively, were set 
equal to zero in the boundary regions of the finite-ele- 
ment model that extend 7/16 in. in from both loaded ends 
of a panel, as illustrated in Fig. 1 . The compression load 
was introduced into the panel by applying a uniform end- 
displacement A to one end of the model w hile holding the 
other end of the panel fixed on the boundary as illustrated 
in Fig. 1: that is, u(0,y) = A and u(L.y) = 0. The simply 
supported boundaries on the unloaded edges (knife-edge 
supports) w-ere simulated by setting the out-of-plane dis- 
placement w equal to zero on a line 5/16 in. in front each 
unloaded edge of a panel. The finite element mesh w as 
defined by using user-written subroutines that are com- 


patible with the STAGS finite-element code. 9 These 
user-written subroutines facilitated the generation of 
models with various cutout sizes and mesh densities, and 
provided a convenient means for assessing the conver- 
gence of a given finite-element model. Convergence 
were done for all finite-element models and a typical 
converged finite-element mesh is shown in Fig. 1. 

Geometrically perfect and imperfect panels were 
studied in the present investigation. Nominal geometry, 
laminate thickness and lamina material properties were 
used for geometrically perfect models. Geometrically 
imperfect models included measured initial geometric 
imperfections, average measured panel thickness and 
thickness-adjusted material properties. The measured 
initial geometric imperfections were included in the 
model by applying a corresponding perturbation to the 
initial geometry of the model by using user-written sub- 
routines. The lamina material properties were adjusted 
for each model using the rule of mixtures and the appro- 
priate average measured panel thickness. 

The predicted nonlinear response of the panels was 
determined with the STAGS finite-element code by 
using a combination of an arc-length path following 
method to model the quasi-static portions of the 
response and a transient analysis method to model 
unstable portions of the response associated with mode 
shape changes that occur during a buckling process. A 
typical finite-element model contained approximately 
10,000 degrees of freedom and a typical nonlinear anal- 
ysis required about 3,000 CPU seconds on an SGI Ori- 
gins 2000 work station. 

Experiments 

A typical curved panel specimen with a centrally 
located circular cutout is shown in Fig. 2 and its corre- 
sponding geometry is shown in Fig. 1. The panels that 
were tested in the study had the same nominal dimen- 
sions and laminate properties as the models described in 
the previous section. The cutouts were machined in the 
panels using a diamond-tip machine tool. The loaded 
edges of the panels were machined flat and parallel in or- 
der to provide uniform end loading during the tests. 

Initial geometric imperfections were measured on 
one side of the panel in order to characterize the actual 
mid-surface geometry of each specimen. The surface 
measurements were taken over a uniform grid that con- 
sisted of 0.125-in. increments along the panel x-axis and 
0.125-in. increments along the panel circumference. The 
maximum magnitude of the measured imperfections, for 
all panels, w as on the order of 10% of the nominal panel 
thickness. In addition, panel thickness measurements 
were taken at several points on the panel to establish the 
average thickness of each panel. 

The panels were loaded in compression by apply- 
ing a uniform end-shortening to the top curved edge of 
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each specimen. The loaded edges were clamped and the 
unloaded edges were supported by knife-edges to simu- 
late simply supported edges. Electrical resistance strain 
gages were used to measure strains, and direct current 
differential transformers (DCDT's) were used to mea- 
sure axial end-shortening and selected out-of-plane dis- 
placements of a panel. In addition, a shadow moire 
technique was used to monitor out-of-plane displace- 
ments and the observed moire fringe patterns were re- 
corded by using video and still photography. All 
specimens were loaded to failure. 

Results and Discussion 

Numerically predicted and experimentally mea- 
sured results for the compression-loaded 15- and 60-in- 
radius panels with a circular cutout are presented in this 
section. The predicted results that are presented subse- 
quently were obtained from finite-element models that 
include the appropriate measured initial geometric im- 
perfections, unless it is specified otherwise. First, results 
are presented for the 15-in-radius panels and then similar 
results are presented for 60-in-radius panels. Finally, re- 
sults that show the effects of elastic circumferential load- 
ed-edge restraint on the prebuckling and buckling 
response of a geometrically perfect 60-in-radius panel 
with no cutout are presented to explain some anomalies 
that were encountered in the test-analysis correlation. 
The values of axial load P and various displacement 
quantities, presented herein, are normalized wdth respect 
to the corresponding linear bifurcation buckling load of 
a panel without a cutout P° cr , that was obtained from a fi- 
nite-element analysis, and the nominal panel thickness t 
= 0.12 in., respectively, unless it is specified otherwise. 
In addition, values of axial stress resultant N x are nor- 
malized with respect to the corresponding linear bifurca- 
tion buckling stress resultant N x ° cr which is equal to P° cr / 
W. The values of P° cr that were obtained for the 15- and 
60-in-radius panels are 66.49 kips and 17.44 kips, re- 
spectively. 

15-in-Radius Panels 

Typical predicted and measured load-shortening 
response curves for the 15-in-radius panels with various 
size cutouts are shown in Fig. 3. Two groups of seven 
curves are show n in Fig. 3 that correspond to values of 
0 < d/W < 0.6 . The solid curves correspond to the 
measured experimental results and the dashed curves 
correspond to the predicted results. Buckling points are 
marked by filled squares and circles for the measured 
and predicted response curves, respectively. An “x” in- 
dicates the failure of a test specimen. 

The measured load-shortening response curves for 
the 1 5-in-radius panels show good correlation with the 
corresponding predicted responses, for the full range of 


cutout sizes. More specifically, the maximum difference 
between predicted and measured values of initial preb- 
uckling stiffness is 3% and the maximum difference be- 
tween predicted and measured buckling load values is 
59c , w'hich occurs for the panel with d/W = 0.1. For the 
panels with d/W = 0.1 and 0.2. however, the unstable 
buckling events resulted in catastrophic failure of the 
panels and thus resulted in no residual postbuckling 
strength. These two catastrophic failures appear to be 
caused by material compression-type failures that initiat- 
ed near high strain concentrations on the free edge of the 
cutout and then propagated almost instantaneously to the 
outer edge of the panel during the failure event. The pan- 
els with values of d/W = 0, 0.3, 0.4, and 0.5 have various 
levels of residual postbuckling strength and can thus sus- 
tain additional load before failure occurs. For the panel 
with d/W = 0, the panel failed near the intersection be- 
tween the loaded-boundary fixture and the knife-edge 
support because of a high strain concentration in that re- 
gion. For values of d/W = 0.3, 0.4, 0.5, and 0.6, the panel 
failures appear to be caused by an interlaminar shear fail- 
ure near the cutout that is caused by large out-of-plane 
bending deformations that occurred near the cutout. 

The results in Fig. 3 indicate that a cutout has a sig- 
nificant effect on the prebuckling, buckling and post- 
buckling responses of the compression-loaded quasi- 
isotropic curved panel considered herein. In particular, 
the predicted and measured results indicate that a signif- 
icant increase in cutout size results in a significant de- 
crease in the overall effective prebuckling stiffness of the 
panel and a substantial increase in the degree of nonlin- 
earity of the prebuckling load-shortening response curve, 
just like the aluminum panels of Ref. 4. The measured 
and predicted results also indicate that for values of d/W 
ranging from 0 to 0.3, the panels have well defined buck- 
ling points (marked by a filled circle and square symbols 
in the figure) and show an unstable buckling event that 
exhibits a significant reduction in axial load and a dis- 
tinct change in the panel displacement field. However, 
for the panels with values of d/W = 0. 1 and 0.2, the buck- 
ling points are marked with an “x” to indicate that the 
buckling event also coincided with the failure of the pan- 
el. In addition, as the value of d/W increases from 0.1 to 
0.3, the buckling load decreases monotonically. Howev- 
er, the buckling load for a panel with d/W = 0 is about 59 c 
less than that for the panel with d/W = 0.1. For values of 
d/W = 0.4 through 0.6, the results show that the response 
no longer exhibits an unstable buckling event. For the 
panels with d/W = 0.4 and 0.5, the measured results indi- 
cate that each of the panels reach a corresponding limit 
point that is followed by relatively benign decrease in 
load before additional load can be applied to the panels. 
This benign limit-point behavior is accompanied by a 
gradual change in panel displacements instead of an in- 
tense, dynamic change in the displacement field like that 
exhibited by panels with small values of d/W. For the 
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panel with d/W = 0.6. the response is characterized by a 
monotonically increasing load-shortening response 
curve that does not exhibit a buckling event. For panels 
with these cutout sizes, large out-of-plane deformations 
develop near the cutout at the onset of loading and their 
amplitudes grow monotonically in a stable manner 
throughout the loading. The benign behavior that is ex- 
hibited by the panels with 0.4 < d/W < 0.6 is unlike the 
behavior shown for flat plates with the same values of d / 
W (see Ref. 8). In addition, the corresponding linear bi- 
furcation buckling loads no longer have any physical sig- 
nificance since the panels no longer exhibit a buckling 
response nor do they exhibit a sudden change in effective 
axial stiffness associated with a buckling event. 

Numerically predicted and experimentally mea- 
sured effective initial prebuckling stiffnesses and buck- 
ling loads for the 15-in-radius panels are summarized in 
Figs. 4 and 5, respectively. The effective initial prebuck- 
ling stiffnesses are defined as the initial slope of the load- 
shortening curves in Fig. 3. Three groups of results are 
presented in Figs. 4 and 5. The solid curves correspond 
to predicted results from geometrically imperfect mod- 
els, the dashed lines correspond to results from geomet- 
rically perfect models, and the filled squares correspond 
to experimentally measured results. 

The results in Figs. 4 and 5 show good agreement 
between analysis and experiment and indicate a reduc- 
tion in the prebuckling stiffness of about 54% as the val- 
ue of d/W increases from 0 to 0.6. Similarly, the results 
indicate a 53% reduction in the buckling load as d/W in- 
creases from 0.1 to 0.5. Moreover, the results indicate, 
for the most part, that the inclusion of the measured im- 
perfections in the numerical models have essentially no 
effect on the predicted initial prebuckling stiffness for 
0 < d/W < 0.6 and the value of the predicted buckling 
load for 0. 1 < d/W < 0.4 . Predicted results for the panel 
with d/W 7 = 0, however, indicate about a 14% reduction 
in the predicted buckling load because of the initial geo- 
metric imperfections. Moreover, the predicted results 
for d/W 7 = 0 and 0.1, that include measured initial geo- 
metric imperfections, show the same trend as the mea- 
sured results; that is. the panel with d/W = 0.1 has a 
higher buckling load than the panel with d/W = 0. The 
predicted and measured results show about a 4.8% and 
5.3% increase, respectively, in the buckling load as d/W 
increases from 0 to 0. 1 . These results at first glance sug- 
gest that the buckling load of the panel w ith d/W = 0 may 
be more sensitive to initial geometric imperfections than 
the panel w ith a value of d/W = 0.1. However, examina- 
tion of the measured initial geometric imperfections for 
these two panels showed that the panel with d/W = 0 had 
a maximum imperfection amplitude that w'as nearly 
twice that of the panel with d/W = 0.1. 

Measured and predicted prebuckling out-of-plane 
displacement 8 at a point on the edge of the cutout arc 


presented in Fig. 6 where 8 = w(L/2. W/2 + d/2). Two 
groups of seven curves are shown in Fig. 6 that corre- 
spond to 0 < d/W 7 < 0.6 . The solid curves correspond 
to the measured experimental results and the dashed 
curves correspond to the predicted results. These results 
indicate that, as the size of the cutout increases, the non- 
linear prebuckling displacements increase dramatically. 
This trend is consistent with the increase in the nonlin- 
earity of the prebuckling load-shortening response 
curves w'ith increasing cutout size that is shown in Fig. 3. 
The predicted results for panels with larger values of d/ 
W agree well with measured results, however, correla- 
tion between predicted and measured results for panels 
with small values of d/W' is poor. 

Predicted out-of-plane displacement contours and 
observed moire fringe patterns for 15-in-radius panels 
with d/W = 0, 0.2, and 0.4 are shown in Figs. 7, 8, and 9, 
respectively. Because of the curvature of the panels, the 
observed moire fringe patterns that are presented are 
somewhat distorted, however, some common features of 
the predicted displacement contours and the observed 
fringe patterns can be identified. The dashed contour 
lines in the predicted contours represent inward displace- 
ments and the solid lines represent outward displace- 
ments. The density of the contour lines is indicative of 
the severity of the displacement gradients in the panel. 

Moire fringe patterns for a panel with d/W = 0 at 
the onset of buckling and initial postbuckling points are 
shown in Figs. 7a and 7b, respectively. The fringe pat- 
tern at the onset of buckling indicates slight inward dis- 
placements in the upper right corner of the panel. Upon 
buckling, the panel snaps dynamically into a mode shape 
that consists of one half-wave along the panel generator 
and two half-waves across its circumference. The pre- 
dicted postbuckling displacement contours, shown in 
Fig. 7d, correlate well qualitatively with the observed 
fringe patterns. It is difficult, however, to determine if 
the predicted displacement contours at the onset of buck- 
ling, shown in Fig 7c, correlate well with the observed 
fringe pattern because of the distortion effects of the pan- 
el curvature on the pattern and the lower density of the 
fringes. 

Observed fringe patterns for a panel with d/W = 0.2 
at the onset of buckling and initial postbuckling points 
are shown in Figs. 8a and 8b, respectively. Predicted re- 
sults for the buckling and initial postbuckling displace- 
ment contours are shown in Figs. 8c and 8d, respectively. 
The predicted buckling displacement contours agree 
well with the observed fringe pattern in Fig. 8a. These 
results indicate that, at buckling, the panel exhibits rela- 
tively large out-of-plane bending displacements near the 
edge of the cutout. Buckling of this panel resulted in cat- 
astrophic failure of the panel (hence P/P° C1 .= 0), thus the 
fringe pattern is not associated with a postbuckling mode 
shape. It is. however, qualitatively similar to the initial 
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postbuckling displacement contour shown in Fig. 8d. 

The observ ed buckling and initial postbuckling dis- 
placement fringe patterns for a panel with d/W = 0.4 are 
shown in Figs. 9a and 9b. respectively. The observed 
displacements at buckling consist of large outward bend- 
ing deformations near the cutout edge and several inward 
buckles away from the cutout. As loading continues in 
the postbuckling range, the displacement fringe pattern 
in Fig. 9a exhibits a gradual change into the fringe pat- 
tern shown in Fig. 9b that consists of one large inward 
buckle to the left of the cutout w hile retaining a fringe 
pattern on the right side of the panel that is similar to the 
corresponding region in Fig 9a. This observed change in 
the displacement field happens in a stable, benign man- 
ner. The predicted results in Figs. 9c and 9d indicate 
good qualitative agreement between the observed fringe 
patterns and the predicted displacement contours. 

The predicted fringe patterns that are shown in 
Figs. 7-9 exhibit essentially either reflective symmetry 
or antisymmetry about the arc x = L/2 or the line y = W/ 
2. Similar contours for the corresponding geometrically 
perfect panels indicate that the small amount of deviation 
from these symmetries are caused by the small, asym- 
metric initial geometric imperfections and not by flexur- 
al anisotropy. This fact is consistent w-ith the small 
values of the anisotropy parameters (g = d = 0.04), for the 
laminate considered herein, that are presented in Ref. 10. 

Predicted stress resultant distributions on the end (x 
= 0) of a 1 5-in-radius geometrically imperfect panel with 
three different cutout sizes and for various load levels are 
shown in Figs. 10, 1 1 and 12. The stress resultant N x is 
normalized w ith respect to the linear-bifurcation buck- 
ling stress resultant of a corresponding panel without a 
cutout N x ° cr , and the circumferential coordinate y is nor- 
malized w'ith respect to the panel arc-width W such that 
0 < y/W < 1.0 . Three curves are shown in each of the 
figures that correspond to predicted stress resultant dis- 
tributions for panels three different values of d/W. The 
solid curves correspond to results for d/W = 0, the dashed 
curves correspond to results for d/W = 0.2 and the dash- 
dotted curves correspond to results for d/W = 0.4. 

Predicted prebuckling stress resultant distributions 
for panels with dAV = 0, 0.2. and 0.4 at a value of P/P° cr 
= 0.5 are presented in Fig. 10. These results indicate that, 
at this load level, the stress distribution for a panel with 
d/W 7 = 0 is. mostly, uniform and symmetric about the line 
y/W = 0.5. The sawtooth shape at the far edges of the dis- 
tribution is caused by the simulated knife-edge supports 
and is an artifice of the localized bending of a panel near 
the knife-edge supports. As the value of d/W increases 
the load is redistributed towards the edge of the panel in 
a symmetric manner. The predicted stress resultant dis- 
tributions at the onset of buckling for panels with d/W = 
0. 0.2. and 0.4 are shown in Fig. I i. In this figure, the 


onset of buckling corresponds to P/P° cr = 0.77. 0.57. and 
0.40 for panels with d/W = 0. 0.2. and 0.4, respectively. 
These results and the results in Fig. 10 indicate that the 
stress resultant distribution for the panels with d/W 7 = 0 
and 0.2 remains symmetric up to buckling with only an 
overall increase in the magnitude of the stress resultant. 
However, for the panel with d/W 7 = 0.4, the stress result- 
ant distribution changes considerably w hen the buckling 
load is approached, w hich is seen by comparing the cor- 
responding results in Figs. 10 and 11. More specifically, 
the central region of the panel edge becomes practically 
unloaded and the stress resultant distribution become 
asymmetric as the buckling load is approached. Unload- 
ing of the central region of the panel is caused by increas- 
ingly large prebuckling out-of-plane bending 
deformations near the cutout (as shown in Fig. 9c) that 
diminish the axial stiffness of the central region of the 
panel which cause the significant redistribution of the 
load towards the simply supported edges of the panel. 
This load path shift is consistent with the reduction in the 
effective stiffness of the panels that is shown in Fig. 3 as 
loading approaches the buckling load. The loss of sym- 
metry of the stress resultant distribution is a result of an 
asymmetric growth in the nonlinear prebuckling defor- 
mations in the panel that is caused by asymmetry in the 
initial geometric imperfection shape. 

Predicted initial postbuckling stress resultant distri- 
butions for panels with d/W = 0, 0.2, and 0.4, that corre- 
spond to values of P/P° cr = 0.53, 0.43 and 0.39, 
respectively, are presented in Fig. 12. The stress distri- 
butions for the panels with d/W = 0 and 0.2 are signifi- 
cantly different from the corresponding stress 
distribution at the onset of buckling (see Fig. 1 1) because 
of the substantial asymmetry in the initial postbuckling 
displacement fields (see Figs. 7-9) and the redistribution 
in load that accompanies buckling. For these panels, the 
stress distributions exhibit minimum and maximum 
stress values associated with large inward and outward 
buckles in the panel, respectively, as shown in Fig. 7d 
and 8d. The initial postbuckling stress resultant distribu- 
tion for the panel with d/W = 0.4 is similar to the stress 
resultant distribution at the onset of buckling (see Fig. 
11), however, the minimum stress resultant value de- 
creases somew'hat upon the formation of a large buckle 
on one side of the cutout, as seen in Fig. 9d. In addition, 
it is interesting to note that the minimum stress values for 
the panels w ith d/W 7 = 0.2 and 0.4 have reversed sign. 

60-in-Radius Panels 

Numerically predicted and experimentally mea- 
sured load-shortening response curves for the 60-in-radi- 
us panels are presented in Fig. 13. Two groups of seven 
curves are shown in Fig. 13 that correspond to values of 
0 < d/W <0.6. The solid curves correspond to mea- 
sured experimental results and the dashed curves corre- 
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spond to the predicted results. Buckling points are 
marked by filled squares and circles for the measured 
and predicted response curves, respectively. An “x" in- 
dicates the failure of a test specimen. 

The results in Fig. 13 indicate that the measured and 
predicted trends associated with the effects of cutout 
size on the load-shortening response of a 60-in-radius 
panel correlate well for the larger values of d/W. How- 
ever, there is a large discrepancy in the results for the 
smaller values of d/W, especially for d/W = 0. In partic- 
ular, the measured buckling load for d/W = 0 is almost 
30% greater than the predicted buckling load. In con- 
trast, the measured buckling loads for d/W = 0.1 and 0.2 
are greater than the predicted buckling loads by 7% and 
2%, respectively. 

The measured and predicted results in Fig. 
Vindicate the panels with values of d/W = 0 to 0.2 ex- 
hibit a significant amount of prebuckling nonlinearity 
and have well defined buckling points that are marked by 
the filled symbols in the figure. Moreover, the prebuck- 
ling nonlinearity in the load-shortening response curves 
is much more pronounced for the 60-in-radius panels 
than for the 15-in-radius panels. An unstable postbuck- 
ling response occurs beyond the buckling point for each 
panel that is characterized by a reduction in axial load 
and a distinct dynamic change in the panel displacement 
field. In addition, as the value of d/W increases, the 
buckling load decreases monotonically. Moreover, the 
panels sustain additional load that is substantially larger 
than the initial postbuckling load until catastrophic fail- 
ure occurs, which is marked by an “x” in the figure. Un- 
like flat plates, the results indicate that the panels with d / 
W = 0.3 to 0.6 no longer exhibit a buckling event; rather, 
the response is characterized by a monotonically increas- 
ing load-shortening response curve up to failure. The 
predicted and measured results indicate that the panels 
exhibit a transition from a buckling response to a mono- 
tonically increasing stable response between d/W = 0.2 
and 0.3, unlike the results for the 15-in-radius panels 
which exhibit the transition between d/W = 0.4 and 0.6. 
These results indicate that curvature has a significant ef- 
fect on the character of the nonlinear compression re- 
sponse of panels with a central circular cutout. 

Another important point to note is that the predicted 
and measured buckling loads all have values that are 
much greater than the predicted linear bifurcation buck- 
ling load for the 60-in-radius panel without a cutout. For 
example, the panel with d/W = 0 has a buckling load 75% 
greater than the corresponding linear bifurcation buck- 
ling load. The high buckling load values are caused by 
large prebuckling deformations that develop in the panel 
and act to retard the onset of buckling. These nonlinear 
prebuckling deformations are neglected in a linear bifur- 
cation buckling analysis. Moreover, the substantial non- 
linear prebuckling response that is exhibited by the panel 
with d/W = 0 is contrary to known behavioral character- 


istics of compression-loaded shallow curved panels. 
These large differences in the buckling load and corre- 
sponding linear bifurcation buckling load indicate that a 
linear bifurcation analysis of a curved panel may not rep- 
resent accurately the actual buckling behavior of shallow 
curved panels in some cases. 

Two types of failure were observed for the 60-in- 
radius panels. Panels with d/W = 0 to 0.3 exhibited cat- 
astrophic compression failures near the test fixture on the 
loaded ends and displayed little evidence of material fail- 
ures near the cutout. Panels with d/W = 0.4 to 0.6 exhib- 
ited interlaminar shear failures near the free-edge of the 
cutout that was caused by extensive out-of-plane bend- 
ing deformations in the region. 

Predicted and measured effective initial prebuck- 
ling stiffnesses and buckling load values for the 60-in-ra- 
dius panels are presented in Figs. 14 and 15, respectively. 
The effective initial prebuckling stiffnesses in this figure 
are defined as the initial slope of the load-shortening 
curves in Fig. 13. Three groups of results are presented 
in Figs. 14 and 15. The solid curves correspond to pre- 
dicted results from geometrically imperfect models, the 
dashed lines correspond to results from geometrically 
perfect models, and the filled squares correspond to the 
experimentally measured results. 

The results in Fig. 14 show good correlation be- 
tween the predicted and measured trends and indicate a 
reduction in initial prebuckling stiffness of about 48% as 
d/W increased to 0.6. The measured imperfections are 
showm to have no influence on the effective prebuckling 
stiffness of the panels. The results in Fig. 15 show good 
correlation between predicted and measured values of 
the buckling loads for d/W = 0. 1 and 0.2, and poorer cor- 
relation for d/W = 0. Overall, all the results indicate a re- 
duction in buckling load as d/W increases from 0 to 0.2. 
Unlike the results for the effective initial prebuckling 
stiffness, the addition of measured imperfections for d/W 
= 0 results in a 5.6% reduction in the buckling load with 
respect to the corresponding buckling load for the geo- 
metrically perfect panel, while the imperfections for d/W 
= 0.1 and 0.2 result in a 1% and 0.5% increase in the 
buckling load, respectively. This trend, and the fact that 
the amplitudes of the measured initial geometric imper- 
fections for these panels are nearly the same, indicates 
that the buckling load of a panel with a small value of d/ 
W is more sensitive to initial geometric imperfections 
than the panels with the larger values of d/W. Moreover, 
the panels with d/W > 0.2 do not exhibit a buckling event 
and, as a result, do not exhibit imperfection sensitivity. 

A comparison between numerically predicted and 
experimentally measured out-of-plane displacements 5 
at a point on the edge of the cutout is presented in Fig. 16 
where 8 ~ w(L/2, W/2 + d/2). Two groups of seven 
curves arc show'll in Fig. 16 that correspond to 0 < d/W 
< 0.6 . The solid curves correspond to the measured ex- 
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perimental results and the dashed curves correspond to 
the predicted results. 

Like the load-shortening results in Fig. 13, the re- 
sults in Fig. 16 indicate that, as the cutout size increases, 
the nonlinear prebuckling displacements increase dra- 
matically. In addition, it is interesting to note that the 
magnitudes of the predicted and measured prebuckling 
out-of-plane displacements are significantly larger for 
the 60-in-radius panels than for the 15-in-radius panels 
(compare Figs. 3 and 13). This trend is consistent with 
the greater degree of nonlinearity of the prebuckling 
load-shortening response curves that is exhibited by the 
60-in-radius panels, compared to the 15-in-radius panels. 
For values of d/W = 0.5 and 0.6, the predicted and mea- 
sured results correlate well, but as the value of d/W de- 
creases, the correlation becomes increasingly worse. 

Predicted and observed out-of-plane displacement 
contours for d/W = 0, 0.2 and 0.4 are presented in Fig. 
17, 18, and 19, respectively. The dashed contour lines in 
the predicted contours represent inward displacements 
and the solid lines represent outward displacements. The 
density of the contour lines is indicative of the severity 
of the displacement gradients in the panel. 

The observed moire fringe pattern for d/W = 0 at 
the onset of buckling is shown in Fig. 17a. The fringe 
pattern is characterized by an hour-glass-shaped outward 
radial displacement pattern. Upon buckling, the panel 
snaps through into a postbuckling mode shape that con- 
sists of a single, large centrally located inw'ard buckle, as 
shown in Fig 17b. The predicted postbuckling displace- 
ment contours in Fig 17d agree well with the observed 
fringe pattern in Fig. 17b, however, the predicted dis- 
placement contours, just prior to buckling, in Fig 17c do 
not agree well with the observed fringe pattern in Fig 17a. 
The observed moire fringe pattern of a panel with d/W = 
0.2 at the onset of buckling is shown in Fig. 18a. The 
fringe pattern indicates outward out-of-plane bending of 
the cutout edges and inward radial displacements be- 
tween the cutout and the loaded edges of the panel. The 
initial postbuckling-displacement fringe pattern is shown 
in Fig. 18b and consists of an inward buckle pattern that 
is similar to the pattern that is shown in Fig. 17b for d/W 
= 0. The predicted displacement contours that corre- 
spond to the onset of buckling and initial postbuckling 
are shown in Figs. 18c and 18d, respectively, and agree 
well with the observed fringe patterns in Figs. 18a and 
1 8b. The observed moire fringe pattern for P/P° cr = 0.63 
and d/W = 0.4 is presented in Fig. 19a and is character- 
ized by large outward out-of-plane bending displace- 
ments at the cutout edge. The displacement field shown 
in Fig. 19b indicates that the bending deformations near 
the cutout have become increasingly severe as the load- 
ing increases. The predicted results in Fig 19c and 19d 
correlate well with the observed fringe patterns in Fig 
19a and 19b. respectively. 

A comparison of predicted stress resultant distribu- 


tions on the loaded end x = 0 of the 60-in-radius panels 
with d/W = 0, 0.2, and 0.4 at various load levels are pre- 
sented in Figs. 20, 21, and 22. The stress resultant N x is 
normalized with respect to the linear bifurcation stress 
resultant of a panel without a cutout N x ° cr , and the cir- 
cumferential coordinate y is normalized with respect to 
the panel arc-width W such that 0 < y/W <1.0. Three 
curves are shown in each of the figures that correspond 
to predicted stress resultant distributions for panels with 
d/W = 0, 0.2 and 0.4. The solid curves correspond to re- 
sults for d/W = 0, the dashed curves correspond to results 
for d/W = 0.2 and the dash-dotted curves correspond to 
results for d/W = 0.4. 

The predicted prebuckling stress resultant distribu- 
tions for a load level of P/P° cr = 0.5 are shown in Fig. 20. 
These results indicate that, for d/W = 0, the center of the 
panel and the supported edges have a slightly greater 
stress resultant magnitudes than the rest of the panel. As 
the value of d/W increases, the center of the panel be- 
comes unloaded and the load path is shifted towards the 
simply supported edges. 

Predicted stress resultant distributions for panels d / 
W = 0 and 0.2 at the onset of buckling, that correspond 
to values ofP/P° cr = 1.35, and 1.23, respectively, and for 
d/W = 0.4 at P/P° cr = 0.63 (the load level at w'hich the 
first distinct displacement pattern appears) are presented 
in Fig 21. The results in Figs. 20 and 21 indicate that, as 
the loading is increased in the prebuckling range, all the 
stress resultant distributions remain symmetric about the 
line y/W = 0.5. Moreover, the stress resultant distribu- 
tions for panels with d/W = 0.2 and 0.4 remain similar in 
shape while increasing in overall magnitude. For d/W = 
0, however, an increasingly larger portion of the load is 
redistributed towards the center of the panel. This re- 
sponse is caused by an increase in curvature of the panel 
that occurs during loading, as indicated in Fig. 17c, 
which causes a geometric stiffening of the panel and a 
shift of the load path into the stiffer central portion of the 
panel. As the cutout size increases, stiffness is lost in the 
central region of the panel which impedes the load path 
from being centrally located. 

Predicted initial postbuckling stress resultant distri- 
butions for d/W = 0 and 0.2, that correspond to values of 
P/P° cr = 0.73. and 0.74, respectively, and for d/W = 0.4 
at P/P° cr = 0.63 (the load level at which the second dis- 
tinct displacement pattern appears) are presented in Fig 
22. In all cases, extensive deformations in the central 
portion of each panel have diminished significantly the 
overall load-carrying capacity and have caused most of 
the load to be redistributed toward the simply supported 
edges of the panels. In particular, the center of the panels 
with d/W = 0 and 0.2, have unloaded completely. 

The predicted prebuckling, buckling and poslbuck- 
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ling stress resultant distributions for the 60-in-radius, 
geometrically imperfect panels shown in Figs. 20-22 are 
symmetric throughout the load history of the panels. In 
contrast, the predicted stress resultant distributions for 
the 15-in-radius panels (Figs. 10-12) are symmetric at a 

load level of P/P° cr = 0.5 but become asymmetric after 
buckling. These differences in the stress resultant distri- 
butions of the 15- and 60-in-radius panels are a conse- 
quence of the panel curvature. More specifically, the 60- 
in-radius panels exhibit symmetric initial postbuckling 
deformations w hile, in contrast, the 15-in-radius panels 
exhibit asymmetric initial postbuckling deformations. 

Anomalous Behavior 

For the most part, correlation between the numeri- 
cally predicted results and the experimentally measured 
results presented herein is good. However, significant 
discrepancies and unexpected behavioral characteristics 
were identified for some of the panels and warranted fur- 
ther study. Most notably, the 60-in-radius panels exhibit 
a significant amount of nonlinearity in the prebuckling 
portion of the load-shortening response curve that is 
shown in Fig. 13. In addition, the 60-in-radius panels ex- 
hibit buckling load values that are significantly greater in 
magnitude than the corresponding predicted linear bifur- 
cation buckling load for a panel without a cutout. These 
response characteristics are contrary to the known be- 
havior of compression-loaded curved panels such as that 
given in Refs. 1 1 through 13. In all cases, the measured 
buckling loads for 15- and 60-in-radius panels are larger 
in magnitude than the corresponding predicted buckling 
load values. The largest discrepancy between measured 
and predicted buckling load values occurs for the 60-in- 
radius panel with d/W = 0 in which case the predicted 
buckling load value is approximately 30% less than the 
measured value. Furthermore, the predicted and ob- 
served buckling deformations for the 15- and 60-in-radi- 
us panels with d/W - 0 do not agree well. The observed 
discrepancies are generally more pronounced for panels 
with relatively small values of d/W, and they tend to di- 
minish as the value of d/W increases. 

As a first step towards identifying the source of the 
significant prebuckling nonlinearity that is exhibited by 
the 60-in-radius panels, the effects of the circumferential 
boundary conditions v(0, y) = v(L, y) = 0 and N xy (0, y) 
= N xy (L, y) = 0 on the response were examined. The nu- 
merical results for the latter traction-free boundary con- 
dition indicated that allowing the loaded edges to expand 
circumferentially in an unrestricted manner eliminated 
the prebuckling nonlinearity of the response. Bectiu.se of 
this result, the result that the panel w ith v(0, y) = v(L. y) 
= 0 underpredicts the degree of prebuckling nonlinearity 
seen in the lest results, and because the test fixture actu- 
ally applies some restraint to the circumferential dis- 
placements of the loaded edges of a panel, another edge 


condition was considered. This edge condition is the 
case in which the circumferential displacement of the 
loaded edges are elastically restrained. This type of 
boundary condition was used as a "first-order" engineer- 
ing approximation to the actual edge restraint that is 
caused by the test fixture. Elastic circumferential edge 
restraints on the loaded ends of a panel are simulated in 
the finite-element model by applying an externally-gen- 
erated linear circumferential stiffness contribution to the 
nodes in the boundary regions that extend 7/16 in. in 
from both loaded ends of the panel, as illustrated in Fig. 
1. The magnitude of the stiffness component w'as varied 
from 0 Ib/in. to 1.0x10 s lb/in. in order to simulate cir- 
cumferential restraints that range from free expansion to 
fixed conditions. As a point of reference, the circumfer- 
ential membrane stiffness of the panels is given by A 22 = 

9.22xl0 5 lb/in. Results that show some of the effects of 
elastic circumferential edge restraints on the compres- 
sion response of the 60-in-radius panel with d/W = 0 are 
presented in the following section. 

Effects of Elastic Circumferential Edge Restraints 

Load-shortening response curves for a 60-in-radius 
geometrically perfect panel with various values of elastic 
circumferential stiffness K v are presented in Fig. 23 and 
illustrate the overall effects of an elastic circumferential 
edge restraint on the compression response. Five curves 
are shown in Fig. 23 that correspond to values of 0 < K v 
< °° where K v = 0 is equivalent to N xy = 0, and K v = <» 
is equivalent to v = 0 on the loaded edges of a panel. The 
solid curve corresponds to measured experimental re- 
sults and the dashed curves correspond to the predicted 
results. Buckling points are marked by filled squares and 
circles for the measured and predicted response curves, 
respectively. The values of the axial load P and the end- 
shortening A are normalized by the corresponding bifur- 
cation buckling load P° cr of a geometrically perfect panel 
with v(0, y) = v(L, y) = 0, and the nominal panel thick- 
ness t, respectively. 

The results in Fig. 23 indicate that a linear-elastic 
circumferential edge restraint has a significant effect on 
the response of a 60-in-radius compression-loaded 
curved panel. For a value of K v = 0 (N xy (0, y) = N XV (L, 
y) = 0), the prebuckling load-shortening response curve 
is linear and the panel exhibits a buckling point at P/P° cr 
= 0.87. Fora value of K v = l.OxlO 5 lb/in., the load-short- 
ening behavior is significantly different. For this case, 
the panel exhibits a significant amount of prebuckling 
nonlinearity and has a relatively high buckling loud of P/ 

P tr = 1 -93. This increase in the nonlinearity of the pre- 
buckling load-shortening response is caused by an in- 
crease in the out-of-plane prebuckling deformations of 
the panel which result in a corresponding decrease in the 
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instantaneous effective axial stiffness of the panel, but 
w ith an overall stiffening effect, with respect to the buck- 
ling resistance. The results show that the out-of-plane 
prebuckling deformations retard the onset of buckling 
and results in about a 30% increase in the buckling load 
with respect to the predicted buckling load for the panel 
with v(0. v) = v(L, y) = 0 on the loaded boundary. For 
values of K v = l.OxlO 7 Ib/in. and °o , the panel exhibits 
less nonlinearity in its load-shortening response curves 
and the buckling loads decrease to P/P° cr = 1.56 and 
1 .44, respectively, as compared to the results for the pan- 
el with K v = l.OxlO 5 lb/in. It is important to note that the 
magnitude of the buckling load and the severity of the 
nonlinear prebuckling deformations are not bounded by 
the results for K v = 0 and K v = <», which is contrary to 
intuition. In addition, the measured buckling load lies 
between the predicted buckling loads for the panels with 
K v = l.OxlO 5 lb/in. and 1.0x10' lb/in. 

The effects of elastic circumferential edge re- 
straints on the transverse center displacement 5 are illus- 
trated in Fig. 24 where 8 = w(L/2, W/2). Five curves 
are shown in this figure that correspond to values of 0 < 
K v < °° . The solid curves correspond to measured ex- 
perimental results and the dashed curves correspond to 
the predicted results. Buckling points are marked by 
filled squares and circles for the measured and predicted 
response curves, respectively. 

For a value of K v = 0, the prebuckling center dis- 
placement is small and inward. Near the buckling load, 
and in the postbuckling range of loading, the magnitude 
of the inw'ard displacement grows rapidly to many times 
the thickness of the panel. For values of K v = 1.0x10'' lb/ 

in., l.OxlO 7 lb/in. and °° , the center displacement re- 
sponse is characterized by a monotonically increasing 
outw'ard displacement. This result is consistent with the 
increase in the degree of nonlinearity of the prebuckling 
load-shortening response. The results also indicate that 
the experimentally measured center displacement re- 
sponse is initially bounded by the results for = 

l.OxlO 5 Ib/in. and K v = l.OxlO 7 lb/in., but as the loading 
increases in the prebuckling range, the measured results 
gradually move out of the bounds of the predicted re- 
sults. 

The effects of elastic circumferential edge re- 
straints on the predicted circumferential (v) displace- 
ment of one corner of a 60-in-radius panel is illustrated 
in Fig. 25. where v 0 = v(0, 0). Three curves are show n in 

Fig. 25 that correspond to values of 0 < K v < °° . Buck- 
ling points are marked by filled squares. The values of 
circumferential displacement v 0 are normalized by the 
panel thickness t. The vertical dashed line corresponds 
to v 0 = O and is included in this figure for clarity. 


For a value of K v = 0, the corner of the panel ex- 
pands continuously in the negative y- circumferential di- 
rection throughout the compression response. For a 

value of K v = l.OxlO 7 lb/in., the corner of the panel ex- 
pands circumferentially throughout the compression re- 
sponse. but to a much lesser degree than that seen for the 

panel with K v = 0. For a value of K v . = l.OxlO 5 Ib/in., 
however, the corner initially expands circumferentially, 
but then, upon further loading, begins to contract circum- 
ferentially. This reversal in the circumferential displace- 
ment takes place at a value of P/P° cr = 0.94 w'hich 
corresponds closely with the load level at which the load- 
shortening response curve transitions from a linear to a 
nonlinear shape (see Fig. 23) and the panel exhibits a 
rapid growth in the magnitude of transverse center dis- 
placement, as indicated in Fig. 24. 

Predicted displacement contours for panels at the 
onset of buckling with various values of K v are presented 
in Fig. 26. The dashed contour lines in Figs. 26(b) 
through 26(e) represent inward displacements and the 
solid lines represent outward displacements. The density 
of the contour lines is indicative of the severity of the dis- 
placement gradients in the panel. In addition, the buck- 
ling load values are included in the figure for each of the 
panels. 

These results in Fig. 26 indicate that the displace- 
ment response changes significantly as the circumferen- 
tial boundary stiffnesses K v increases. For a value of K v 
= 0, the panel displacements are characterized by a large 
centrally located inward deformation pattern, as shown 
in Fig. 26b. For values of K v = l.OxlO 7 Ib/in. and °° , the 
predicted panel displacement are characterized by a large 
centrally located outward deformation pattern as shown 
in Figs. 26d and 26e, respectively. For a value of K v . = 

l.OxlO 5 lb/in., however, the predicted displacement con- 
tours in Fig 26c are characterized by an hour-glass pat- 
tern similar to the observed moire fringe pattern shown 
in Fig. 26a. 

A summary of the predicted effects of elastic cir- 
cumferential edge restraints on the magnitude of the 
buckling load of a 60-in-radius, geometrically perfect 
panel is presented in Fig 27. Predicted buckling load val- 
ues are marked in the figure by filled square symbols and 
connected with solid line segments. 

The predicted results in Fig. 27 can be placed into 
three distinct groups. The first group of results are for 

4 

panels with 0< K v <3.0.vl() . These panels exhibit a 
linear prebuckling load-shortening response similar to 
the predicted response for K v = 0 that is shown in Fig 23. 
The buckling load values for these panels are all below 
the predicted linear bifurcation buckling load. In addi- 
tion, these panels exhibit transverse center displacement 
responses, circumferential corner displacement respons- 
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es and displacement contours similar to the predicted re- 
sponse for K v = 0 (N xy = 0) shown in Fig 24. 25 and 26b. 
respectively. The second group of results are for panels 

w ith 4.0.v 1 O'* < K v < 3.0a 10^. The panels in this group 
exhibit a significant nonlinear prebuckling response sim- 
ilar to the predicted response discussed previously herein 
for a panel with K v =1.0xl0 5 lb/in. The prebuckling 
load-shortening response for this group of panels exhib- 
its a significant amount of nonlinearity, and the buckling 
load values can exceed twice the predicted linear bifur- 
cation buckling load value. In addition, the character of 
the displacement response is similar to that predicted for 
the K v =1.0xl0 5 lb/in. case that is illustrated in Figs. 24, 
25, and 26c. The final group of results are for panels 

with 4.0a 10^ < K v < °° . The panels associated with this 
group exhibit response characteristics similar to those 
predicted for the fully clamped, K v = °° (v = 0) case, that 
is illustrated in Figs. 23, 24, 25 and 26e. 

Concluding Remarks 

Results from a numerical and experimental study of 
compression-loaded quasi-isotropic curved panels with a 
centrally located circular cutout have been presented. 
The numerical results were obtained by using a geomet- 
rically nonlinear finite-element analysis code. Some of 
the effects of cutout size, panel curvature, initial geomet- 
ric imperfections, and elastic circumferential edge re- 
straints on the prebuckling, buckling and postbuckling 
the response have been described. These results indicate 
that a cutout can have a significant effect on the prebuck- 
ling, buckling and postbuckling response of a compres- 
sion-loaded panel. For example, the results indicate that 
an increase in the cutout size results in a decrease in the 
effective prebuckling stiffness of the panel, as expected, 
but causes an increase in the nonlinearity of the load- 
shortening response curve. For relatively small values of 
the cutout-diameter-to-panel-width ratio d/W, the panels 
exhibit a distinct unstable buckling event. For relatively 
large values of d/W, the panel response no longer exhib- 
its an unstable buckling event; rather, the panel response 
is characterized by a monotonically increasing load- 
shortening cune. The results also show that the value of 
d/W that marks the transition between these two types of 
nonlinear behavior depends significantly on the panel 
curvature. In particular, panels w ith a 15-in-radii exhibit 
the transition at a value of d/W 0.4 < d/W < 0.6, and the 
60-in-radius panels exhibit the transition between d/W 
0.2 < d/W < 0.3. Results have been presented that also 
show that the buckling loads for panels with small values 
ol d/W 7 are more sensitive to initial geometric imperfec- 
tions than panels with large values of d/W. This trend is 
consistent w ith the idea that, as the amount of material in 
the center ol the panel decreases and the prebuckling de- 


formations near the cutout dominate the response, the 
panels become less sensitive to initial geometric imper- 
fections. When the load-shortening curve is monotoni- 
cally increasing and there is no buckling event, there also 
is no imperfection sensitivity. 

Correlation between the numerically predicted and 
experimentally measured results w<as good in most cases. 
However, significant discrepancies and unexpected be- 
havioral characteristics were obtained for certain cases. 
For example, the shallow 60-in-radius panels w'ith small 
values of d/W exhibit a significant amount of nonlinear- 
ity in the prebuckling portion of the load-shortening re- 
sponse curve and exhibit buckling load values that are 
significantly greater than the predicted linear bifurcation 
buckling load for a 60-in-radius geometrically perfect 
panel without a cutout. This behavior is contrary to pre- 
viously known behavioral characteristics of compres- 
sion-loaded shallow' curved panels. Predicted and 
observed buckling deformations for 15- and 60-in-radius 
panels with d/W = 0 do not correlate well. Furthermore, 
predicted buckling load values for 15- and 60-in-radius 
panels with small values of d/W are consistently under- 
predicted when compared to the experimentally mea- 
sured buckling load values. The largest discrepancy 
between measured and predicted buckling load occurs 
for a 60-in-radius panel with d/W = 0, in which case the 
predicted buckling load value is 30% less than the mea- 
sured value. However, the observed discrepancy be- 
tween the measured and predicted buckling loads 
decreases as the value of d/W increases. 

In an effort to explain some of the behavioral issues 
stated above, results from a limited numerical parametric 
study of the effects of elastic circumferential edge re- 
straints on the compression response of a 60-in-radius 
geometrically perfect panel have been presented. These 
results indicate that an elastic circumferential edge re- 
straint on the loaded edges of a curved panel has a signif- 
icant effect on the behavioral characteristics. Overall, 
the study has shown that the predicted results can be 
placed in to three groups that exhibit distinct behavioral 
characteristics depending on the magnitude of the elastic 
circumferential restraint. Panels in the first group with a 
relatively low' magnitude, approximately 3% of the pan- 
els circumferential membrane stiffness exhibit a lin- 
ear load-shortening response that contains a distinct 
buckling point. The magnitude of the buckling load for 
these panels is alw ays less than the predicted linear bifur- 
cation buckling load and the prebuckling deformations 
are negligible. Panels in the second group have an inter- 
mediate magnitude of the elastic circumferential re- 
straint, 4% of A 22 to 3.25A22, exhibit a significant 
amount of nonlinearity in the prebuckling response and 
can exhibit buckling loads that are more than tw'ice the 
predicted linear bifurcation buckling load. Panels in the 
third group are for relatively high magnitudes of circum- 
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ferential restraint, greater than 3.25A22- exhibit some 
nonlinearity but not nearly as pronounced as the nonlin- 
ear behavior exhibited by the second group of panels. 
The predicted results for the panels with an intermediate 
magnitude of the elastic circumferential restraint more 
closely represent the experimentally measured results 
and indicate that some amount of restricted circumferen- 
tial motion at the loaded boundary is present in the test 
fixture during the experiments. 
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Fig, 1 Typical model geometry and boundary conditions 
(dashed lines mark the rows on which the boundary condi- 
tions are applied). 
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Fig. 2 Typical test specimen and apparatus. 
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Fig. 3 Numerically predicted and experimentally mea- Fig. 5 Numerically predicted and experimentally mea- 
sured load-shortening response curves for 15-in-radius sured buckling load for 15-in-radius panels with a cutout, 

panels with a cutout 



d/W 

Fig. 4 Numerically predicted and experimentally mea- 
sured effective prebuckling stiffness for 15-in-radius pan- 
els with a cutout. 
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Fig. 6 Numerically predicted and experimentally mea- 
sured out-of-plane center displacements for 15-in-radius 
panels with a cutout. 



(a) Moire fringe pattern at the onset ol buckling, (b) Initial postbuckling moire fringe pattern, 

P/P° cr = .77 P/P° cr = .53 

Fig. 7 Observed and predicted out-of-plane displacements for a 15-in-radius panel with d/W = 0. 
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(c) Predicted displacements at the onset of 
buckling, P /P° cr = .57 


td) Predicted initial postbucklinp displacements, 
P/P 1 ' = .43 


Fi}j. 8 Observed and predicted out-of-plane displacements for a 15-in-radins panel with d/W = 0.2. 


13 

American Institute of Aeronautics and Astronautics 














y/W 


Fig. 12 Initial postbuckling stress resultant distribution on 
the loaded end of a 15-in-radius panel with a cutout. 



d/W 

Fig. 14 Numerically predicted and experimentally mea- 
sured effective initial prebuckling stiffness for 60-in-radius 
panels with a cutout. 



Fig. 13 Numerically predicted and experimentally mea- Fig. 15 Numerically predicted and experimentally mea- 
sured load-shortening response curves for 60-in-radius sured buckling loads for 60-in-radius panels with a cutout, 

panels with a cutout. 



2 r Measured response 


Predicted respons 



lig. 16 Numerically predicted and experimentally measured out-of-plane displacements for 15-in-radius panels with a cutout. 
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(c) Predicted displacements at the onset of (d) Predicted initial postbuckling displacements, 

buckling, P/P° cr = 1.35 P/P° cr = 73 

Fig. 17 Observed and predicted out-of-plane displacements for a 60-in-radius panel with d/W = 0. 



(a) Moire fringe pattern at the onset of buckling, (b) Initial postbuckling moire fringe pattern, 

P/P" cr = 1.26 P/P° cr = .75 

Fig. 18 Observed and predicted out-of-plane displacements for a 60-in-radius panel with d/W = 0.2. 
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(c) Predicted displacements at the onset of (d) Predicted initial postbuckling displacements, 

buckling, P/P° cr = 1.23 P/P° cr = .74 


Fig. 18 Continued. 



(a) Moire fringe pattern at P/P° cr = .63 (b) Moire fringe pattern at P/P° cr = 1.6 



(U Predicted displacements at P/P“ t . r - .63 (d ) Predicted displacements at P/P" 1T = 1.6 

lig. 19 Observed and predicted out-of-plane displacements for a 60-in-radius panel with d/VV = (1.4. 
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Fig. 20 Stress resultant distributions on the loaded end of a 
60-in-radius panel with a cutout and P/P° cr = 0.5. 


Fig. 21 Stress resultant distributions on the loaded end of a 
60-in-radius panel with a cutout at the onset of buckling. 
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Fig. 22 Initial postbuckling stress resultant distributions 
on the loaded end of 60-in-radius panclwith a cutout. 



Fig. 23 Effects of circumferential boundary stiffness K v on 
the load-shortening response of a 60-in-radius panel with 



Fig. 24 Effects of circumferential boundary stiffness K v on 
the out-of-plane center displacement of a 60-in-radius 
panel with d/W = 0. 



Fig. 25 Effects of circumferential boundary stiffness K v on 
the circumferential displacement of a corner of a 60-in-in- 
radius panel with d/W = 0. 
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Fig. 27 Effects of circumferential boundary stiffness K v on the buckling load of a 60-in-radius panel with d/W = 0 
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